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We further pursue an approach to the sign problem of quantum chromodynamics at nonzero 
chemical potential, in which configurations of the lattice path integral are gathered into subsets. 
In the subset construction we multiply each temporal link by center elements independently and 
in a first step neglect the gauge action. The positivity of the subset weights - shown for 0+1 
dimensions in an earlier study - extends to larger lattices: for two sites in the temporal direction 
and arbitrary spatial extent we give a proof of the positivity by decomposing the subset weight in 
positive summands. From numerical evidence we conjecture that the positivity persists on larger 
lattices and that the gauge action can be reintroduced through mild reweighting. First results on 
the quark number obtained with this method in two dimensions are shown as well. 

PACS numbers: 12.38.Gc 


I. INTRODUCTION 

The sign problem in quantum chromodynamics (QCD) 
- the non-positivity of the quark determinant at nonzero 
chemical potential - hampers numerical simulations of 
QCD thermodynamics based on importance sampling, 
see e.g. [T]. This problem, generically caused by a com¬ 
plex action in the partition function, also occurs in other 
physical systems ElE]. The sign problem causes large 
cancellations between the contributions of the configura¬ 
tions corresponding to the fundamental degrees of free¬ 
dom of the ensemble, and suggests that other degrees of 
freedom may be more effective to render the path integral 
in such regimes. 

The subset method consists of gathering configurations 
of the ensemble into subsets using some definite rule. 
The weight of such subsets is given by the sum of the 
individual weights. An ‘early application of subsets’ (to 
summations) is the computation of finite arithmetic se¬ 
ries by first pairing up numbers positioned symmetrically 
around the center of the sequence and then summing up 
the pair sums (this goes back to the Indian mathemati¬ 
cian Aryabhata [1] and is also attributed to the young 
C.F. Gauss). In lattice QCD and related theories the 
aim is to find subsets with positive weights. The remain¬ 
ing path integral over all subsets can then be treated by 
statistical methods like importance sampling. 

The subset method has been first developed on a ran¬ 
dom matrix model of QCD: subsets with positive weights 
could be found, hence, solving the sign problem 
Later, this subset construction was understood in terms 
of imaginary chemical potentials and canonical partition 
functions [7] . Canonical partition functions vanish, when 
the corresponding configurations average out in the inte¬ 
gral. The subset method removes such contributions al¬ 
ready on the level of the integrand, rendering the weight 
positive. 
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Although gathering conhgurations in subsets and 
adding up their weights is a very general and always exact 
option to compute partition functions, there is generically 
no reason that the resulting subset weight be positive. 
In this context it is very helpful, if the generation of the 
configurations in the subset is guided by some symmetry 
principle or physical insight to perform cancellations ex¬ 
plicitly using more physical degrees of freedom. For the 
random matrix model this symmetry is a .Zjv subgroup 
of the U(I)-symmetry of the integration measure. For 
lattice QCD, multiplying any of the group-valued links 
by another group element is possibl^ i.e. compatible 
with the Haar measure. A natural choice, however, is 
to make use of the center symmetry as a subgroup 
of the gauge symmetry SU(3), on the temporal lattice 
links0 Mathematically speaking we divide out the cen¬ 
ter from these links. This leads to temporal mesonic and 
baryonic hoppings. As another motivation we take the 
fact that gauge groups without the center, like SU(A^) 
in the adjoint representation or G 2 , do not cause a sign 
problem [SHII]. 

In our previous work with T. Wettig [5] center sub¬ 
sets have been shown to yield positive weights in one 
of the simplest QCD toy models, QCD in O+I dimen¬ 
sions, where only the Dirac determinant is present (and 
where analytic solutions are available [mia). Since the 
configurations in this system are fully characterized by 
one group element, the Polyakov loop, the subset method 
consists of a single .Za-multiplication. This again removes 
canonical determinants that average out in the partition 
function, and only multiples of the baryon chemical po¬ 
tential survive. 

In the present work we extend the subset method to 
higher dimensions, using (unrooted) staggered fermions 


^ Multiplying gauge links by generic group elements has been uti¬ 
lized on top of the center subsets in [S] to treat many flavors. 

^ Changing a single link and thereby changing the configurations’ 
weights, as we do here, should not be confused with a gauge 
transformation that always changes several links and preserves 
all gauge invariants. 
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and in a first step neglecting the gauge action. We 
demonstrate that the method yields positive weights, 
sums of Dirac determinants, when each temporal link is 
multiplied by an independent center element. We prove 
the positivity of these center subset weights analytically 
for lattices with two sites in the temporal direction and 
arbitrary spatial extent and present numerical evidence 
that positivity persists for larger lattices. We show first 
numerical evidence that the gauge action can be rein¬ 
troduced and treated through reweighting. As a first 
measurement using this method we present the quark 
number for the massless case on two-dimensional lattices 
of various sizes. 

As it stands the subset method requires a cost growing 
exponentially with the volume. In the positivity proof, 
however, positive subterms can be identified. This could 
be the starting point of a refined subset method. In a 
companion paper m we give a diagrammatic represen¬ 
tation of the subset weights shedding more light on the 
terms causing the sign problem and how the latter is 
solved by the subset method. 

This paper is organized as follows: in the next section 
we introduce the general idea of subsets and show that 
they improve the reweighting factor in general. Sec. |III| 
contains the definition of subsets for lattice QCD, its ba¬ 
sic properties and building blocks. In Sec. |IV| we prove 
the positivity of the subset weights for lattices with two 
sites in the temporal direction by expressing the fermion 
action in terms of Grassmannians, and after that for the 
massless 2x2 lattice using Dirac matrix language. We 
also comment on the possibility to apply subsets on spa¬ 
tial links. In Sec. |V] we present our numerical results, 
which support the conjecture that this positivity also 
holds for larger lattices, and show data for the quark 
number density in two dimensions. We also give first re¬ 
sults for the subset reweighting factor in the presence of 
a gauge action. Finally, we summarize and give several 
technical results in appendix. 


where |G| is the cardinality of G (the number of elements 
in G and thus in for all x). We divide by this number 
to avoid a |G|-fold overcounting of the configurations in 
the integral. 

In the best case, the new weights are positive for 

all subsets fix- It means that having performed part of 
the integration (the summation over y in fix) explicitly 
and deterministically, the remaining integral (over fix) 
may be subject to importance sampling methods. 

Even if not becoming strictly positive, the integrand 
always ‘comes closer to positivity’ in the sense of an im¬ 
proved sign quenched reweighting factor. Besides the 
invariance of the measure we only need the Cauchy- 
Schwarz inequality to show that 



For the reweighting factors r and their variances we 
obtain the inequalities 


^subsets — 


> 


Jdx\a{flx)\ Jdx\f{x)\ 


— ^sign quenched 5 


/\2 A 2 

^r, subsets — ^r, sign quenched’ 


(4) 


where the second inequality follows from the first one 
and the fact that = 1 — (in this case) [6]. 
This improvement is in some sense expected, since sum¬ 
ming/integrating over the whole ensemble yields the par¬ 
tition function, which shall be positive. 

Note that this reduction of the sign problem comes at 
the expense of an increase of the computational effort by 
a factor of |G|. 


II. IDEA OF SUBSETS AND IMPROVEMENT 
OF THE REWEIGHTING FACTOR 

We consider a general integral Z = Jdy{x) f{x) whose 
integrand f{x) is real, but not necessarily positive. The 
idea is to collect several configurations x and add up their 
weights f{x) to a new weight. To formalize this we as¬ 
sume that the integration measure is invariant under the 
action of a discrete group G, f dy(x) f{gx) = J dy,(x) f{x) 
for all g € G. For the integration over gauge groups in 
lattice QCD, dg,{U) is the Haar measure, which obeys 
this invariance. The subsets fl generated by such a group 
and their weights a read, 

^x = {gx\gGG}, (1) 


HI. SUBSETS FOR LATTICE QCD 
A. Definitions 

A lattice QCD configuration is given by the tuple of 
temporal and spatial SU(3) links, which we denote as 
U = [Uo,Us]- The partition function 

Z = Jd[Uo,Us] det Di[Uo,Us]) (5) 

is a path integral with Haar measure d[UQ,Us] = 
Ha; dg,H{Uo{x)) dfj.H(Ui^x)) . Herein we work in the 
strong coupling limit and neglect the gauge action (in 
Sec.|V]we reintroduce the gauge action through reweight¬ 
ing). We consider the staggered Dirac operator for one 
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quark flavor with mass m and chemical potential /i, 

Di[Uo,U,];x\y) 

d-1 

+ ~ (y)^x-ly] +2TOl34,y, (6) 

i=l 

where we have set the lattice spacing to unity, a = 1, 
and neglected a factor 2 on the left hand side, since in 
the determinant of D this only gives an irrelevant con¬ 
stant factor. We work on an W x x • ■ • x Nd-i lattice 
with even iV^ and antiperiodic boundary conditions in 
the temporal direction (represented by the minus signs 
on the last time slice). The temporal extension is the 
inverse temperature, Nt = 1/T, while the factors 
are the lattice implementation of the chemical poten¬ 
tial m- We use the standard notation 0 and i for unit 
steps in the temporal and spatial directions. For sim¬ 
plicity we have chosen the staggered signs as 770 = 1 and 
rji = (_i) 2 :o-i-a;i+...+a:i_i ^ ^gj-Q ghemical potential the 
determinant of the staggered Dirac operator is positiv^ 
but at nonzero real y the Dirac operator is no longer an- 
tihermitian, as D{yL)^ = and its determinant is 

no longer necessarily positive, inducing the sign problem. 

In the subset method for QCD we generate subsets 
by gathering a number of different configurations of the 
ensemble summing up their individual weights to the 
subset weight. The members of a subset are generated 
by multiplying the temporal links Uq{x) with the cen¬ 
ter Z 3 = {1, z, z*} where z = exp(27rj/3). The generated 
links C/q, zUq and z*Uo remain in the configuration space, 
the group SU(3). The invariance of the Haar measure 
under group multiplications ensures that the integration 
measure is the same for the three configurations, so that 
adding the determinants gives the correct weight of the 
subset (in the strong coupling limit). 

B. 0+1 dimensions 

In 0+1 dimensions the configurations are fully charac¬ 
terized by the (untraced) Polyakov line P (which can be 
shown using a gauge transformation). In Ref. [8] it was 
shown that the subset Up = {P,zP,z*P} for one quark 
flavor has a positive weight for any P G SU(3), 

cr(Dp) = i [det D{P) + det D{zP) + det D{z*P)\ 

= - 3A + A| tr Pp + 2 cosh(3/r/T) > 0 , (7) 


® The massless Dirac operator anticommutes with (the residual 
chiral symmetry), such that its eigenvalues come in diA pairs. 
Antihermiticity at /i = 0 yields eigenvalues on the imaginary 
axis and the positivity of the determinant follows. The mass 
term only shifts the real part of the eigenvalues keeping this 
positivity. 


where A = 2 cosh (Vj arsinh (^ ^)) > 2, thus solving 
the sign problem. In the massless case this simplifies to 
A = 2 and 

cr(flp) = 2 + 2|trP|2+ 2cosh(37i/T) >0 (m = 0). 

( 8 ) 

Note that only terms with baryon chemical potential 
/if, = 3/i appear. 

C. Higher dimensions 

In higher dimensions each temporal link can be multi¬ 
plied by an independent center element. In other words, 
the subsets are generated by the direct product of all 
local Zzix), 

^[Uopis] = {[g^oMsWg G G}, 
y 

G = (g) Z^{xi ), V = NtN ^... Nd-i , (9) 

i=l 

or a subgroup thereof. A group element g = 
^g27rifei/3^..., 2} acts on a con¬ 

figuration Uq = (C/o(2:i), ■ • •, Uo{xv)) through 

gUo = ..., /^Uo(xv)) ■ ( 10 ) 

After introducing subsets the partition function © can 
be written as 

Z = Jd[Uo,Us]a{n[u„,u.]), (11) 

with subset weights 

1 ^ " 

f^i^[Uo,Us]) = wHE detDi[glto,l(s]), (12) 

2^1 ki—0 

where we just add the determinants because of the in¬ 
variance of the Haar measure under group multiplica¬ 
tion. The cardinality of G, |G| = 3^^, is exponential in 
the number V of lattice points. 

D. Collective subsets 

Before discussing the full subsets ([^, let us first con¬ 
sider collective subsets, which contain three configura¬ 
tions. The subset elements are constructed by syn¬ 
chronously rotating the links on one time slice with the 
same Z 3 element exp(27r7fc/3), k G {0,1,2}, and leav¬ 
ing all other links untouched. For the Dirac operator 
^ the operation of a collective Z 3 rotation can also be 
interpreted as adding an imaginary chemical potential 
p/T -G p/T + 27rzfc/3 while keeping the original links 
Uo{x) unchanged, since the chemical potential can be 
introduced equivalently through factors 
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on one time slice. The partition functions with these 
shifted chemical potentials are identical, which is the 
Roberge-Weiss periodicity m- Nevertheless, for a given 
configuration the integrands in the three partition func¬ 
tions differ and after adding these up the contributions 
of the canonical determinants with nonzero triality are 
removed, as we now show. 

The fugacity expansion of the partition function, Z = 
with canonical partition functions Zq, has a 
corresponding expansion at the level of the determinants 
and subset weights, 


deiD{[UoMs])=Y.Dq 

Q 


aW/T 


(13) 


where the sum ranges over the number of spatial lattice 
points, q = —314, ■ • ■, 314, with 14 = Ni... Nd-i, and 
we have omitted the link arguments of Dq and aq. 

As the Z^ rotations can be shifted into an imaginary 
chemical potential, the collective subsets yield for every 
q (see also [7]) 


collective «/i/T _ ^ 




^qiti,IT+2Trik/3) 


(14) 


fc=0 


and the well-known formula for the sum over powers of 
roots of unity reduces q to multiples of three, 

DqC / - E ^q{27Tik/3) _ ^q^ijT 


/c=0 


^collective 


= D, S, 


q Oq^Sb 


(15) 

with integer baryon number b in the range — 14 ,..., I 4 . 
Therefore, the canonical weights for collective subsets are 
given by the canonical determinants with zero triality 
and vanish otherwise. It is well-known that the QCD 
partition function satisfies the same property, i.e. only 
receives triality zero contributions, and can therefore be 
expanded in the baryon chemical potential fXb = 3/r. One 
of the essential mechanisms of subsets is that they ensure 
this reduction already at the level of the path integrand. 
Although the collective subsets and the ensuing global 
reduction to triality zero terms turned out to solve the 
sign problem in O-fl dimensions [S], we will see in Sec.[V| 
that in higher dimensions they attenuate the sign prob¬ 
lem but do not suffice to solve it. Further note that the 
collective subsets preserve all plaquette values, such that 
this subset method could be applied directly in full QCD 
after introducing the gauge action. 

For the full subsets Q the cancellation of the collective 
subsets is also achieved (since the collective subsets are a 
subgroup of the full group), but the zero triality weights 
will be modified further. Full subsets enforce zero triality 
for each temporal link, i.e. even for chemical potentials 
that would be defined locally, again via /r(a;) —?► yt{x) + 
2TTik{x)/3. This statement will be made more precise in 
the next section, e.g. Eq. (21 1 . 


E. Subset building blocks 


In this subsection we will demonstrate how the sub¬ 
set weight can be decomposed into local building blocks. 
They are the basics for the positivity proof presented in 
the next section. For that we write the quark determi¬ 
nant in the partition function as an integral over Grass- 
mann fields. 


det D{U) 


■ y 3 

nn dlpa{Xi)d'lpa{Xi) 

J-1 a=l 


(16) 


where S'f(i/', V', ^) = the 

staggered Dirac operator of Eq. §■ As each term in the 
fermion action is bilinear in Grassmannians, the expo¬ 
nential can be factorized in single-link contributions 

w) = exp {'ilj{x)D{U;x\y)-ip{y)) . (17) 

x,y 


Let us focus on the contribution of a specific tempo¬ 
ral link, i.e. Uo{x) and Uq{x) for some fixed x, which 
according to § is 

e^'il}{x)UQ{x)4’{x -f 6) — e~^'i(}{x + 6 )Uq{x)'iP{x) 

= e^u + e~^u , (18) 

which defines the bosonic variables u and u. At the tem¬ 
poral boundary this expression enters the fermionic ac¬ 
tion with a minus sign. We expand the contribution from 
this link to the path integral weight 

^ 1 

exp{ef^u + e-f^u) = V -^(e'^M)"(e-^i 2 )“. (19) 

n\ ml 

n,m—0 


Because V’ and ip represent Grassmann fields with three 
color degrees of freedom, only terms which are at most 
cubic in u and u will contribute. 

Center rotations on that link amount to the changes 
Uo{x) —)■ e'^'^^^^^Uo{x), and thus u —>■ and u 

g-2irifc/3y^ k g { 0 , 1 , 2 }, and the subset sum for a term 
(e^u)”(e“'^u)'" becomes 




gMe2^ife/3y)n(g-Fg-2,rifc/3-)^ 




27r2A;(n—m)/3 


g(n niod 3 = 0, 

0 otherwise. 


( 20 ) 

( 21 ) 


Terms with nonzero triality are removed, while terms 
with zero triality remain unchanged. Therefore, the fol¬ 
lowing subset building block^survive from the local con- 


^ Note, that the symbol a now represents subset sums over Grass¬ 
mann terms which are still subject to Grassmann integration. 
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tribution (19): 


aj; = l + uu+ ^{uu)'^ + ^{uu)^ + + e 

( 22 ) 


We will denote the term, just three forward hoppings, 
by ‘baryonic’, the term, just three backward hoppings, 
by ‘antibaryonic’, the terms with (mw)", n G {1, 2, 3}, by 
‘n-mesonic’ (hoppings) and the identity by ‘empty link’. 
Only the (anti)baryonic terms carry the /i-dependence, 
in the form , with baryon chemical potential 


l^b- 

The subset building blocks are similar to those found in 
the polymer-baryon approach to strong coupling lattice 
QCD [T7l[T8]. In the latter the gauge links are integrated 
out completely, yielding the same hopping structure (in 
all directions) as in Gx- Integrating out the Grassmanni- 
ans as well, enforces site constraints on the polymers and 
baryons. In contrast to this, our hopping terms retain a 
dependence on the background gauge links (actually only 
the 1- and 2-mesons do, see (29)), but are subject to the 
same Grassmann constraints. 


IV. POSITIVITY 

In this section we first present our proof for the positiv¬ 
ity on lattices with Nt = 2 and arbitrary spatial extent. 
Then we give a more detailed discussion of the 2x2 lat¬ 
tice, where some concepts become more explicit. While 
the first subsection uses the framework of Grassmanni- 
ans, the second subsection makes use of the particular 
matrix structure of the Dirac operator (which is of course 
fully equivalent). 


A. Positivity proof for Nt = 2 using Grassmannians 


We first consider the temporal part of the fermionic 
action and gather the two points in the temporal direction 
at each spatial point. 


SF,t = ■ 


(23) 


where in analogy to 


u= ^p{l,x)Uo{l,x)^p{2,x), 
u = -tl!{2,x)U^{l,x)tl}{l,x), 
v= '0(2,f)f7o(2,f)'(/)(l,f), 

V = -1/5(1, x)U^ (2, x)tp(2, x ), 


(24) 


with arguments x on the left hand sides omitted for sim¬ 
plicity. The relative minus sign between u and i; in (23) 
comes from the antiperiodic boundary conditions. 

For the subset sum we expand exp{SF,t) in analogy to 
Sec. |IIIE[ and again the zero triality condition applies to 
{u,u) and {v,v). Using (22) and taking into account the 


additional sign for v and v the contribution of the spatial 
point X to exp(S'F_() is 





r u 

■ 

( 




X U —)■ —V, U 


1-b 

-v), 


liT 


uu - 


(uuY 

2(2 


(25) 


where four of the hopping terms in (22) were collected 


into a product. Because of periodicity, u and v visit the 
same Grassmann variables, although they are connected 
by different gauge links, and the same holds for v and u. 
The Grassmannian antisymmetry restricts (25) to poly¬ 
nomials that are cubic both in the combinations (it, v) 
and {u,v). Therefore, several cross terms from the prod¬ 
uct above vanish and one is left with 


^x — ^(2,a;) — ? 


(26) 


where 


(7b = 




o 






= -3ai 


pS/a 


3! “ 3! ^ 


(27) 


combines baryonic terms with the empty link and 3- 
mesons, and the remainder 


Cm = 1 -I- UM -I- 


(uuY 

2 ! 


{uu —>■ vv) — 1 


X (uu 


[uuY {vvY 

= uu + vv-\—^ -1— ^ —I- uu vv 

uu(vvY (uuYvv 


2 P 


2(2 


(28) 


is a polynomial in uu and vv and thus mesonic. Note that 
various contributions in a have Grassmann vacancies and 
need to combine with spatial hoppings or mass terms to 
achieve Grassmann saturation for all i/’(x). 

Using the antisymmetry of Grassmannians and 
det 17 = 1 one can sirnplify the baryonic terms to be 
independent of the link£[j 


^ = ~ = -Y^(l,d:)Y^(2,x), 
^ = ~ = -Y^( 2 ,x)Y^(Yx), 


(29) 


^ We can write: 

= YYYU'ff = YYYaUa_a.'Va')iYbUw<Pb'){YcU^FVc') 

^ Y-abc^a'b'c'^aa't-^bb'U^ 

= — ^ (3! det 

where itaYbYc = ea6cV’i’/’2V’3 = e-abcY^ with i/;® = ipiip2'4>3, and 
we used the Leibniz formula for the determinant and det U = 1. 
The minus signs comes from the permutation of Grassmann vari¬ 
ables. 
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where = ipi'4’2'4’3- The baryonic product can therefore 
be written as 


CTb = 


1 + 2cosh(3^) 


3! 


1 + 2cosh(3^) 


3! 


= [l — 2 cosh(3^)'0^(l, a:)] 

X [l + 2cosh(3/r)^/i^(2,;?)'0^(l,ai')] . 


(30) 


With (30) and (28) the subset sum (26) for the tem¬ 


poral links at each site x can be written as 

as = [l — 2 cosh(3/r)'(/j^(l, x)4’^{2, a?)] 

• [l — 2 cosh(3^)^^(l, al)Vi^(2, a:)] 


U-U + V- V+ — ■ — + — ■— +UV-UV 


-2 - 2 

UV UV 


2 - -2 

U V U V 


2 ’ 


(31) 


where we have rewritten the terms such that the first 
factor depends on {■0(1,af),'!/'(2,a:)} and the second on 
{^(1,af),-0(2,a;)}, cf. the definitions (24) of u,u,v and 


V (for clarity the multiplication of such factors is rep¬ 
resented by a dot). With appropriate functions ha the 
subset weight ([sT]) can be summarized as 


= T), 0(2, x),Uo{l, x), Uo{2, x)*) 


a—O 

•h„(0(l,f),0(2,a-),t/o(l,a-)M/o(2,f)). (32) 

In each term both factors can be obtained from one 
another by exchanging the fermions and antifermions, 
0^0, and complex conjugating the links. After denot¬ 
ing this exchange operation as 


/ = / 




we can write 


as = '^ ha{x) ■ ha{x). 


(33) 

(34) 


For the spatial part of the fermion action, see ( |17[ ), we 
find, after anticommuting the fermions in the backward 
hopping term. 


Sf,s = X! X! {^a.{x)r]i{x)Ui^ab{.x)'ii}b{x + i) 

X i 

+ il}{x)ar]^{x)Ul^t^{x)-^b{x -f i 
= (w(a;) -h w(a;)) , 


(35) 

(36) 


where we explicitly wrote out the color indices in (35) 


(with implicit summation over repeated indices) to easily 
identify the exchange symmetry (33). The spatial weight 


(without subsets for the spatial links) is then 


= JJexp {w{x)) ■ exp (^w{x)). (37) 


This embodies the antihermiticity of the Dirac operator, 
which is an important ingredient in the positivity of the 
determinant at /i = 0, and hints at the fact that the 
structure (34) will be important to prove the positivity 


of the subset weights. 

The second ingredient necessary to show the positiv¬ 
ity of the subsets is the staggered chirality, i.e. the fact 
that all interactions connect even with odd sites or odd 
with even sites (recall that m = 0). Here even and odd 
lattice sites x^ and x° are those with 775 ( 3 ;) = - 1-1 and 
— 1, respectively, where 775 = (—all 
products ha - ha in (34) and e’" • e'^ in ( [^ one factor 
will only depend on the Grassmann sets { 0 ( 3 :®)} = 
and { 0 ( 3 ;°)} = 47 ° coming with links {U^{x‘^)} = W 
and {Ul{x°)} = U°* and the other only on the com¬ 
plements 47 ®, with links 77®*, U°. The full massless 
subset weight reads 

a = ■ ^“s(0) ■ exp(^(*))> 

X 0^1 X 

(38) 

and after expanding all products we can rewrite a as 
a = ^/A(;F®,^°,77®,77°*) -/^ (4-®, ^°,77®*,77°) , (39) 

A 

where the index A = {as} runs over all combinations of 
as for all spatial points x and oversaturated terms, i.e. 
with too many Grassmannians, automatically cancel. 

Moreover, the Grassmann measure factorizes analo¬ 
gously 

Y[di;{x)di;{x) = Y[Y[ d0(x®)fi0(3:®)d0(x°)d0(3;°) 

X X^ x° 

= n d'4){x'') d0(3;°) • d0(a;®) d{[i{x°) 

X^ x° X^ x° 

= (i4/®d^7° . d4'®d4'° , (40) 

where d'tpdij} = Ila^i'^0a'^0a- The Grassmann integrals 
over a remove all terms with insufficient Grassmann con¬ 
tent and each of the surviving terms is a product of two 
polynomials in the (complex) links only. By renaming 
integration variables it can be seen that these polynomi¬ 
als are still related by complex conjugation of the link 
arguments 

f Y[d^p{x)d^p{x)a=J 2 [ (^0 4'°,77®,77°*) 

■ j d4'®d§°/A (4'®,4'°,77®*,7/°) 

= ^p^(77®,77®*)p^(77®*,77°). (41) 

A 

Since the polynomials p emerge from the hopping terms 
in the Dirac operator, they possess real coefficients only. 
Therefore, polynomials with complex conjugated link ar¬ 
guments turn into complex conjugated polynomials, such 
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that 


W_d'il){x)d'il}{x) a = '^PA {WM°*) [PA {WM°*)]* 

X A 

= > 0 - (42) 


This proves that the massless subset weight is positive 
for = 2 at nonzero chemical potential and arbitrary 
spatial extent. What is more, the subset weight consists 
of various positive subterms (labeled by A). 

The modification of this proof caused by a nonzero 
mass is presented in Appendix For p = 0 the subset 
construction is not needed and one can prove the posi¬ 
tivity by using ( |T^ for all directions. 

In showing the positivity we have made use of the fact 
that for Nt = 2, u and v connect the same sites and thus 
have the same Grassmann content but with opposite p- 
dependence, and that the baryonic factors are indepen¬ 
dent of the connecting links, (29), such that p enters in 


the form (301. This would not hold with nonzero triality 


terms present, as is the case in the determinant formula¬ 
tion without subsets. 

Let us make some further remarks. Firstly, note that 
the temporal antiperiodicity is crucial to ensure the sub¬ 
set positivity. Consider for instance periodic boundary 
conditions. In this case the v- and u-terms in (23), (25) 


and (27) would have plus signs instead of minus signs. 


and all odd powers of v and v in the derivation above 
would have opposite signs. The mesonic discussion is 
left unchanged as v and v always come in pairs, but the 
baryonic product would become 




—3/i 




I -I- 2sinh(3^) — 


1-f 


o 3/x 


g3/i 


1 - 2sinh(3/r) — 


(43) 


Due to the different signs in both factors, this product 
cannot be written as h ■ h and no longer satisfies the 
conjugation symmetry required in the positivity proof 
above. 


Note that the last two terms of the mesonic weight (28) 
can be simplified further as (see Appendix 

uu{vv)'^ = {uu)‘^vv 

= 4 I tr P(x)p^/i^(I, a:)^^(l, x)'ijj^(2, af)^/>^(2, x). (44) 


This and the 4cosh^(3p,)-term from (30), which equals 
2 -I- 2cosh(6/r) = 2 -|- 2cosh(3^/T), have full Grass¬ 
mann content and represent the full weight for the one¬ 
dimensional massless case at fVt = 2, cf. (1^. 


The last remark is slightly more formal. One can 
rewrite the subset contribution a as an exponential of 


an effective subset action, 

E / _ _ {uuY + (uu)^ 

\uu -\- vv — ^ 

X 

— Q'ii>'^{l,x)ip'^{l,x)^^{2,x)'ip^{2,x) (45) 

-I- 2cosh(3^){^^(l,x)'!/'^(2,f) - '(/)^(2,f)'(/)^(l,f)}^ , 

which can easily be checked by an expansion of the ex¬ 
ponential function (which again terminates). In contrast 
to the Dirac action this effective action is not bilinear in 
('i/'. Ip) but also involves higher powers of the Grassmann 
fields. Therefore the Grassmann integral of a can not be 
represented as a determinant. 


B. Alternative proof for the massless 2x2 lattice 


Below we give an alternative positivity proof, which 
only holds for a 2 x 2 lattice in the massless case. The 
salient feature of this proof is that it directly uses the 
determinant formulation. 

For a 2 X 2 lattice the Dirac operator can be written 
as 


(0 

0 

Tf 

5i\ 

0 

0 

52 

'Y'e 

-^2 

-TO 


0 

0 

V-^l 

'Y'o 

^2 

0 

0 / 


where and St are the temporal and spatial hoppings 
on the spatial slice x and time slice t, respectively. These 
3x3 blocks are given by 


r° 

rpO 

-^2 

51 

5 2 


e>^Uii{ 2 l) + , 

- e'^17o(22) + e-'^C/o^(I2)l , 

- e'^C/o(12) + e-^C/o^(22)] , 

-' 1 / 1 ( 11 )-C7j(I2)] , 
Ui{22)-U\{2l), 


(47) 


where the superscripts e and o stand for even-odd and 
odd-even hoppings. Each entry is a sum of two contribu¬ 
tions because neighboring sites on a 2 x 2 lattice can be 
connected in two ways, where one is ‘around the world’. 
The spatial part of the Dirac operator is antihermitian, 
therefore only two independent St occur. At ^ = 0 the 
full Dirac operator is antihermitian and T° = (T®)f. The 
sign difference between and S 2 is due to the staggered 
phase, whereas the signs in reflect the antiperiodic 
boundary conditions. 

Using the determinant formula for block matrices we 
find 


det D = det 




= det(S'S'1')det(l3 


r-Tp 

-T^J 

M")det(l3- 


M°), 


(48) 




























where we defined 


5 = S1S2 , 

^-1- 


-1 


= T^S:^^T^S- 
M° = [S:[^]'<T°[S2^]'<T° . 


(49) 


The first factor of (48) is the determinant of the spatial 


part of the Dirac operator, which is blind to ^ and there¬ 
fore positive. Formally, the positivity follows because a 
matrix product SS^ is always positive-semidefinite. At 
fi = 0 the M-matrices are related as M° = {M^)\ and 
the full determinant is positive for the same reason. 

The last two determinants of 3 x 3 matrices in ( [4^ can 
be expanded as 

det(l 3 - M) 1 - det M - tr M -b 

(50) 

which can, for example, be proven in terms of the eigen¬ 
values of M. Before constructing the full subsets, we 
first construct ‘coarse subsets’ flc containing three con¬ 
figurations by multiplying ?7o(ll) with the three center 
phases and 17o(21) by their complex conjugate. The re¬ 
maining links are left untouched. These rotations form a 
subgroup of the full subset group. They multiply Tf and 
thus by the three center phases and Tf and thus M° 
by the complex conjugate phases, and consequently 

det(l 3 - = 1 - det M® - tr M® 


-b 


4 ^»fc /3 (trM®)^-trM 


re 2 


det(l 3 - = 1 - detM° - tr M° 

4,,fe;3(trM°)^-trM°^ 


-b e 


(51) 


The coarse subset sum of the product (48) can then be 


computed as in (211 and we find 


^ det T» = det(S'S'l') 


(1-det M®)(l-det M°) 


tr M® tr M' 


„ (tr M®)2 - tr M® 2 (tr M®)^ - tr M° 


(52) 


After completing the full subsets the chemical potential 
cancels in the second line of (52), since it can only enter 


if a temporal link or its inverse appears three times, and 
by inspection of (47) and (49) this cannot happen. Since 


this expression is independent of /i, its value can equally 
well be computed at /r = 0 where M° = (M®)1. This 
gives the following full subset contribution for this second 
line 

= I det S'ptTo(^| tr M®|^ -b (tr M®)^ — tr M®2|^^ 

(53) 


where cro( - ■ ■) indicates the subset at /r = 0, and clearly 
cr(//) is positive. Full subsets of the first line of (52) give 
a contribution 


(^(i) = X! det(5'5'''') 1 - det(S' 
l-det(S'-3)^detTfdetT2° 


idetTi®detT2® 


(54) 


where we substituted M®’°, defined in ( |49| . The subset 
sums for temporal hoppings are computed in Appendix 
[C|and substitution of (C6) yields. 


(T(/) = I det 5'|2 -b (det S + det 5^) 4cosh^(3^) 
-b4(2cosh2(3/i) -b |trP(l)|2)(p(l) ^ F(2)) 

= |4cosh^(3^)-b debs'll 
-b8cosh2(3/r)(|trP(l)|2 -b |trP(2)|2) 
-b4|trP(l)|2|trP(2)|2, (55) 

where the Polyakov loops P only depend on the spatial 
argument after tracing. This expression contains mesonic 
(/r-independent) and baryonic (3/r-dependent) terms. 

The full subset weight is simply 




lUo Ms 


— cr(/) -b <7(11) ) 


(56) 


with a(j) and cr(jj'f given in (55) and (53). As all the 
summands are positive the subset weight is positive too. 
It is interesting to note that the first term explicitly com¬ 
bines a /i-independent and a /r-dependent contribution to 
achieve its positivity. A similar principle is at work in the 
first term of (311 in the more general proof given in Sec. 

IIY3 

Finally, let us have a look at the canonical subset 
weights (integrands of the subsets canonical partition 
functions) as introduced in Eq. (13). This means nothing 
but collecting the terms in the subset result according to 
their /^-factors. The terms _ g±6/i/T from 

the term 16cosh^(3/i) = (e^^ -b in (56) and have 

weights iJq=±Q = 1 as expected from (antijbaryon sat¬ 
uration of the lattice. The weights of the next terms 

g±6M ^ g±3M/T g^j.g 


(57) 


<7q=±3 — det S -b (det S)* -b 4 -b 2| tr TiP -b 2| tr P 2 \ 

To compute the first two terms we note that 
det S 2 = det ( Ui{22)Ui{21) -l) det Ul (21) 

1F(2) 1 

= tr W{2) - tr 1T(2)1' = 2i Im tr 1F(2), (58) 

where W{t) is the Wilson loop on time slice t closing 
around the spatial boundaries, and the second equality 
can easily be checked for any SU(3) matrix in term of its 
eigenvalues. Similarly we find for the other time slice 


det S'! = —2ilmtr 1F(1), 


(59) 
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where the minus sign comes from the staggered phase. 
On multiplying we find 


det S' = 4Im tr VF(1) Imtr Vl^(2). 


(60) 


The first two terms of (571 can thus be negative and are 
not necessarily compensated for by the remaining posi¬ 
tive terms in aq=±^. Therefore, the subset positivity does 
not generically hold for the individual canonical subset 
weights. 


C. Note on spatial subsets 


In this work we apply center subsets to temporal links 
only, which we conjecture are the minimal subsets for 
achieving positivity of subset weights. The main motiva¬ 
tion for this choice was that the chemical potential caus¬ 
ing the sign problem only couples to the temporal hop¬ 
pings. Extending the idea and applying center subsets 
on the spatial links as well is another option to compute 
the partition function, which we briefly comment on. 

As the temporal subsets are positive already, the spa¬ 
tial subsets just add positive numbers and thus remain 
positive. Obviously, the cost for such subsets is even big¬ 
ger, namely which is why we have not used spatial 
subsets in practice. 

The usefulness of spatial subsets can be seen at the 
level of the canonical subsets, which are not necessarily 
positive as discussed at the end of the previous section. 
Consider the term (IfiOl) for the 2x2 lattice, which can 


cause a negative contribution to crg=±3 in (57). This term 
is linear in the Wilson loops W{1) and W{2) and will 
disappear after spatial subsetting according to Eq. (21) 
(with (n,m) = (1,0) or (0,1)). In this particular case 
even the canonical subset weight becomes positive upon 
spatial subsetting. 

In future work we plan to use the positive summands 
obtained after subsetting to sample the partition func¬ 
tion with a worm algorithm. In this context the spatial 
subsets could give an additional advantage as they fur¬ 
ther reduce the number of allowed building blocks, while 
still keeping the dependence on the SU(3) links. 


V. NUMERICAL RESULTS 


We implemented the full subset method in a Monte 
Carlo simulation. The subset weights are computed ex¬ 
plicitly by adding the numerically computed determi¬ 
nants for all the configurations belonging to the sub¬ 
set. These positive subset weights are then used to gen¬ 
erate relevant subsets of the partition function using a 
Metropolis importance sampling algorithm. 

As the full subset contains an exponential number of 
configurations we speed up the computation in a number 
of ways: 


• There is some redundancy in the full subset: although 

by definition the full subsets contain configura¬ 

tions, there is a 3^*“^ fold degeneracy in the determi¬ 
nant values so that we effectively only need to consider 
^(Vs-i)Nt+i different configurations per subset. 

This degeneracy occurs when all temporal links in¬ 
side time slice i are rotated by the same Za-factor 
Zj, i = while the Polyakov lines are left 

unchanged, i.e. = ^- These constrained rota¬ 

tions lead to a 3^*“^ degeneracy of each determinant 
value. This equality of determinants can be under¬ 
stood as the latter consist of closed loops and the 
rotations leave all loop values unchanged: for loops 
involving temporal links we either encounter Zi and 
z* if the loop goes backward and forward in time, 
albeit at different spatial points, or we wrap around 
the lattice which also yields unity because of the con¬ 
straint. 

• Rank-6 corrections are used to reduce the numerical 
work in the computation of the determinants when 
stepping from one configuration to the next in the 
subset. 


• The algorithm is efficiently parallelized by evenly dis¬ 
tributing the configurations of each subset over several 
threads. 

More details on the numerical implementation will be 
given in a forthcoming publication. 

In Table [T] and Fig. we present results obtained 
with direct product subsets for QCD in two dimensions 
with massless staggered quarks. We compare the aver¬ 
age reweighting factors for two conventional reweighting 
schemes in the link formulation: 


(a) phase-quenched reweighting, and 

(b) sign-quenched reweighting. 


with the reweighting factors for three different subset 
constructions: 


(c) a collective subset constructed by synchronous Z3 
rotations of all temporal links on one time slice, as 


mentioned in Sec. Ill D 


(d) a direct product of local Z 3 subsets for the temporal 
links of all spatial sites on one time slice (containing 
3'^“ configurations), which we call T-slice subsets, 
and 


(e) the full subsets (§, which is a direct product of local 
Z 3 subsets for the temporal links on all lattice sites. 

For the full subsets data were collected for Ntx2 grids 
with Nt = 2, 4, 6,8,10, for At x 4 grids with Nt = 2, 4, 6 
and for a 2 X 6 and 2x8 grid, all for Nf = 1 and 
m = 0 in the strong-coupling limit. As can be seen 
from the phase-quenched reweighting factor (a) the sign 
problem steadily grows as Nt and is increased. The 
sign quenched reweighting (b) somewhat reduces the sign 
problem, which can be useful for simulations at small 
chemical potential [H]. Whereas a collective Z^ rota¬ 
tion (c) does not bring much improvement in the two- 
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Nt X A, 

2x2 

4x2 

6x2 

8x2 

10 X 2 


a 

phase-quenched 

0.8134(3) 

0.4361(4) 

0.233(2) 

0.130(2) 

0.071(1) 


b 

sign-quenched 

0.9271(2) 

0.6150(5) 

0.355(3) 

0.203(2) 

0.109(2) 


c 

collective 

0.9778(9) 

0.777(4) 

0.500(6) 

0.303(8) 

0.178(3) 


d 

T-slice 

1.0 

0.9896(5) 

0.885(2) 

0.670(5) 

0.436(8) 


e 

full 

1.0 

1.0 

1.0* 

1.0* 

1.0* 



Nt X Ntt 

2x4 

4x4 

6x4 

2x6 

2x8 

a 

phase-quenched 

0.7934(5) 

0.295(1) 

0.0961(9) 

0.7364(6) 

0.6725(7) 

b 

sign-quenched 

0.9197(3) 

0.442(2) 

0.149(1) 


0.8917(4) 

0.8523(5) 

c 

collective 

0.959(1) 

0.557(6) 

0.214(8) 


0.912(3) 

0.867(2) 

d 

T-slice 

1.0 

0.9973(2) 

0.812(3) 


1.0 


1.0 


e 

full 

1.0* 

1.0* 

1.0* 


1.0* 


1.0* 



TABLE I. Reweighting factors for 2d-QCD for Nf = 1 (m = 0) for (a) phase-quenched and (b) sign-quenched reweighting in 
the link-formulation, and for (c) collective subsets, (d) T-slice subsets and (e) full subsets. The columns give the data for NtX2 
grids with Nt = 2, 4,6, 8,10, for a A’t x 4 grid with Nt = 2,4, 6 and for a 2 x 6 and 2x8 grid, all in the strong-coupling limit 
a.t fi — 0.3 with A^mc = 100,000 (* means A^mc = 1, 000). 


Nx 2 468 



FIG. 1. Reweighting factors r for 2d-QCD for Nj = 1 (m = 0) in the strong coupling limit at /r = 0.3, where the sign problem 
is largest, with Nmc = 100,000 (for the larger lattices the full subsets were simulated with Nmc = 1,000, see Table |I|. We 
compare the reweighting factors in the link-formulation for phase-quenched (pq) and sign-quenched (sq) reweighting, and in 
the subset formulation for collective, T-slice and full subsets. Each column shows the different (color coded) reweighting factors 
for a specific Nt x N^. lattice. For Nt = 2 the T-slice data are not visible as they overlap with the full direct product data 
(r = 1.0 for both of them). 


dimensional case, the T-slice subsets (d) substantially 
improves on the sign problem. However, the truly sur¬ 
prising observation is that the full subsets (e) yield subset 
weights that are reaj^and positive in all cases considered. 


The imaginary part of the weights is trivially cancelled by im¬ 
plicitly pairing each configuration with its complex conjugate. 


We have proven this property for Nt = 2 in Sec. |IV[ but 
conjecture that it holds for any lattice size in any dimen¬ 
sion. 


As an application of the subset method we show the 
quark number density as a function of T and ^ (in lattice 
units) in Fig.j^for QCD in 1-1-1 dimensions (for the com¬ 
putation of observables in the subset framework, see [SI 
Eq. (3.4)]). One observes the Silver Blaze phenomenon 
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FIG. 2. Surface plots of the quark number density n ver¬ 
sus temperature T and quark chemical potential (in lattice 
units) for Nx = 2 {Nt = 2 ... 12) and Nx = 6 {Nt = 2 ... 8). 


where the quark number is independent of /i below some 
value fic when T —0. For the larger lattices the full 
subsets (e) were too costly and we used T-slice subsets 
(d) instead. These subsets require additional reweighting 
(away from Nt = 2), but still yield a vast improvement 
over the standard phase quenched reweighting method. 

We also verified the effect of the gauge action on the 
reweighting factors for the 2x6 lattice by switching on /3 
to leave the strong-coupling regime. The subset weights 
have to be modified to take into account the different 
values of the gauge action for the different subset ele¬ 
ments, and the sign problem slowly reappears even for 
the full product subsets. Nevertheless, for /3 = 1, 2, 3,4, 5 
the reweighting factor is 1.0, 1.0, 0.984(7), 0.964(13), and 
0.972(17) respectively, so that the sign problem remains 
very mild, at least for these parameter values. 

As a further test we also looked at the full subsets 
for small lattices in three- and four-dimensional QCD, 
even though the computational cost is huge even for 
small lattices. We observe with great interest that for 
2^, 4 X 2^, and 2'^ lattices the full subsets always give 
positive weights, as was verified on samples of 200 ran¬ 
dom configurations. 


VI. SUMMARY 

We have applied subsets generated by center multipli¬ 
cations on temporal links to QCD at nonzero chemical 
potential and have proven the emerging subset weights 
to be positive for Nt = 2 lattices in the strong cou¬ 
pling limit. We have also presented numerical evidence 
that leads us to conjecture that the positivity persists 
for larger lattices. Moreover, preliminary results show 
that the reweighting induced when reintroducing the 
gauge action is not severe. In a subset measurement of 
the quark number density the typical Silver Blaze phe¬ 
nomenon is clearly visible. 

The number of determinants constituting the full sub¬ 
sets grows exponentially with the number of lattice sites, 
i.e. with inverse temperature and volume, which is a rein¬ 
carnation of the sign problem. This is why our numerical 
studies have been restricted to small lattices so far. Two 
ways out of this situation are conceivable. First, smaller 
than full subsets improve the reweighting factor consid¬ 
erably - for which we gave an analytic argument as well 
as numerical evidence - and these subsets are thus help¬ 
ful for extending the applicability range of reweighting 
methods. 

The second possibility relies on our finding that the full 
subset weight can be decomposed in a sum of positive 
terms. Some of their weights simplify considerably, as 
will also be discussed in a further publication M- More¬ 
over, the Grassmann nature of the fermions constrains 
the combinations of such building blocks. Such con¬ 
strained systems can typically be simulated using worm 
algorithms, and will be the subject of future work. 

Open issues are the conjecture about the subset posi¬ 
tivity for Nt > i lattices, which still needs to be demon¬ 
strated, and the need to investigate how large the sign 
problem becomes on larger lattices when introducing the 
gauge action and leaving the strong coupling regime. 
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Appendix A: Positivity of the massive subsets 


In this section we extend the positivity proof given in 
Sec. II V Al to the massive case. 

The mass insertions, up to three per site. 


exp {2m'tp{x)'ij}{x)) =1-1- 2m^{x)ilj{x) 


( 2 to )^ 

2 ! 


{il)[x)^{x)f 


(2m)^ 

3! 


{'ip{x)il){x)f , (Al) 


break the chiral symmetry, since they connect sites with 
themselves, and thus modify the positivity argument of 
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Sec. IV A which was based on the staggered chirality of 
the Dirac operator. 

For a given subset term let n denote the total num¬ 
ber of mass insertions on even sites. The total number 
of mass insertions at odd sites is n toc0 and thus the 
partition function only contains even powers of m. For 
the argument it is irrelevant how many insertions come 
from any given site, the only thing that matters is that 
the coefficients of all even mass terms in (All have the 


same sign and those from all odd mass terms too. We 
denote the mass insertion locations as x° and x® with 
i,j = 1,... ,n and regroup the fermions as follows, 


Indeed, W = U~^ in SU(3) and the left hand side of (B3) 
is the matrix of cofactors, which equals U~^ transposed, 
up to the determinant of U, which is unity in SU(3). 
Substituting (B3) in (Bl) and (B2| yields 


u^v = = 2trP'^ 


c' c cc' 

U‘‘v = —2(p^ip^Uc'cVcc' = —2tTP(p^tf}^ 


(B4) 


with Polyakov line P = UV. Gathering these results 
yields 

= dtrPtrP^(B5) 


n n 


Z=1 


n 

n 


2=1 


n 

2=1 

n 


n n 

[Y[^ixj)][Y[Hxj) 






n n 


2=1 




(A2) 


which comes with a factor These mass insertions 

multiply expressions similar to (38) for the hopping con¬ 
tributions without changing their structure, i.e. they re¬ 
main products of complex conjugate polynomials after 
Grassmann integration. Because the even powers of m 
have positive coefficients and no additional minus signs 
are picked up by reordering the Grassmann variables (a 
reordering of ip is always accompanied by the same re¬ 
ordering of "(/j) the positivity proof holds as before. 


Appendix B: Mixed-meson contribution 


In order to simplify the mixed-meson hopping 
{uu)‘^(vv) we substitute ip = ip{l,x), ip = ip{2,x), 
U = Uii{l,x) and V = Uo(2,x) in (24), for simplicity, 
and expand 


U^V = (ipUpfi-lpV'^p) = Ipa'P’b’lpcUaa'Ubb'Pa'Pb'Pc'V^c' 

= Ip ip dabc^a'b'c'Uaa'Ubb'^cQf 

U^V = {-pU^1pY{(pVlp) = -ipa'Pb'PcUl^.Ulylpa'Ipb'Ipc'Vcc' 
= -'P^1p^eabc(^a'b'c'Ul^,Ul^,Vcc', (B2) 


where we also reordered the Grassmannians. We will also 
use the following relation for U S SU(3), 

2^abc^a'b'c'Uaa'Ubb' — ^c'c' 


^ For n mass insertions on even sites the remaining {3VI2 — n) 
fermions must be provided by temporal and spatial hop¬ 

pings to saturate the Grassmann integrals. These hoppings also 
contain the same number of 1 ^( 1 °). The same is true for 
and 'ip{x°). In order to saturate the Grassmann integrals on the 
odd sites, n mass insertions are required there as well. 


where the sign vanishes after commuting Grassmann 
variables. 


Appendix C: Temporal hoppings for 2x2 lattices 


Herein we compute the subsets of temporal hoppings 
needed in (54). Their definitions (47) as sums of two 
temporal links can be compactly written as 


o _ 

with P = X 


t° = {x + 1) mod 2 . 


(Cl) 


The determinants are computed in the following way 
detT®’° = (-)"+! det(e^C/o(t"’°,x)t7o(t°’",x) + e-'^) 

P(P’°,x) 

X det x) 

' V-^ 

1 

= (-)“+! -h e-'^ tr P{x) + tr P(x)1' -b , 

(G2) 

where the Polyakov loops P only depend on the spatial 
argument after tracing, and the second equality can be 
derived for any P S SU(3) using its eigenvalues. The 
product of these determinants at the same spatial posi¬ 
tion X, 

det T® det T° = -b 2 tr P{x) 

-be"^'^ [2trP(x)^ + (trP(x))^] -b 2(1-b | trP(x)|^) 

-b [2 tr P{x) + (tr P(x)'^)^] -b 2 tr P{x)^ + , 

(G3) 


is nothing but the one-flavor determinant of massless one¬ 
dimensional QGD for W = 2 [HI Eq. (A.6) with A = 2]. 

What is needed in (54) are total subsets on various 
products of such determinants. Since the latter depend 
on the links only through Polyakov loops, it is sufficient 
to center rotate the P{x), which removes their nonzero 
triality terms. From (C2) we obtain the subset as the 
zero triality projection 


cr (det T®’°) = (-)“+! 2 cosh(3Ai) . (G4) 
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Note that T® and T° depen d on the same P{x) such that 
the subset on the product (C3), 


case in 0+1 dimension, cf. Eq. 
on different x factorize, giving 


Finally, the subsets 


CT (det det r°) = + 2 + 2| tr P{x)\^ + e®'" 

= 4cosh^(3/r) + 2| trP(x)p , (C5) 

is not the product of individual subsets. Again, this ex¬ 
pression agrees with the subset weight for the massless 


cr (det det T^) = a (det T° det T^) = —4 cosh^(3/r), 
a(detri"detT2"detT°detT2°) (C6) 

= (4cosh^(3/r) + 2|trP(l)p) x (P(l) ^ P(2)) . 
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